Plasma SARS-CoV-2 RNA Levels as a Biomarker of Lower Respiratory Tract SARS-CoV-2 Infection in Critically Ill Patients With COVID-19.

Plasma SARS-CoV-2 viral RNA (vRNA) levels are predictive of COVID-19 outcomes in hospitalized patients, but whether plasma vRNA reflects lower respiratory tract (LRT) vRNA levels is unclear. We compared plasma and LRT vRNA levels in serially collected samples from mechanically ventilated patients with COVID-19. LRT and plasma vRNA levels were strongly correlated at first sampling (n = 33, r = 0.83, P < 10-9) and then declined in parallel in available serial samples except in nonsurvivors who exhibited delayed vRNA clearance in LRT samples. Plasma vRNA measurement may offer a practical surrogate of LRT vRNA burden in critically ill patients, especially early after ICU admission.

Plasma SARS-CoV-2 RNA levels as a biomarker of lower respiratory tract SARS-CoV-2 infection in critically ill patients with COVID-19.

Plasma SARS-CoV-2 viral RNA (vRNA) levels are predictive of COVID-19 outcomes in hospitalized patients, but whether plasma vRNA reflects lower respiratory tract (LRT) vRNA levels is unclear. We compared plasma and LRT vRNA levels in simultaneously collected longitudinal samples from mechanically-ventilated patients with COVID-19. LRT and plasma vRNA levels were strongly correlated at first sampling (r=0.83, p<10 -8 ) and then declined in parallel except in non-survivors who exhibited delayed vRNA clearance in LRT samples. Plasma vRNA measurement may offer a practical surrogate of LRT vRNA burden in critically ill patients, especially early in severe disease.

Severe Acute Respiratory Syndrome Coronavirus 2 Viremia Is Associated With Coronavirus Disease 2019 Severity and Predicts Clinical Outcomes.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral RNA (vRNA) is detected in the bloodstream of some patients with coronavirus disease 2019 (COVID-19), but it is not clear whether this RNAemia reflects viremia (ie, virus particles) and how it relates to host immune responses and outcomes. METHODS: SARS-CoV-2 vRNA was quantified in plasma samples from observational cohorts of 51 COVID-19 patients including 9 outpatients, 19 hospitalized (non-intensive care unit [ICU]), and 23 ICU patients. vRNA levels were compared with cross-sectional indices of COVID-19 severity and prospective clinical outcomes. We used multiple imaging methods to visualize virions in plasma. RESULTS: SARS-CoV-2 vRNA was detected in plasma of 100%, 52.6%, and 11.1% of ICU, non-ICU, and outpatients, respectively. Virions were detected in plasma pellets using electron tomography and immunostaining. Plasma vRNA levels were significantly higher in ICU > non-ICU > outpatients (P < .0001); for inpatients, plasma vRNA levels were strongly associated with higher World Health Organization (WHO) score at admission (P = .01), maximum WHO score (P = .002), and discharge disposition (P = .004). A plasma vRNA level >6000 copies/mL was strongly associated with mortality (hazard ratio, 10.7). Levels of vRNA were significantly associated with several inflammatory biomarkers (P < .01) but not with plasma neutralizing antibody titers (P = .8). CONCLUSIONS: Visualization of virus particles in plasma indicates that SARS-CoV-2 RNAemia is due, at least in part, to viremia. The levels of SARS-CoV-2 RNAemia correlate strongly with disease severity, patient outcome, and specific inflammatory biomarkers but not with neutralizing antibody titers.

CpG Methylation Profiles of HIV-1 Pro-Viral DNA in Individuals on ART.

The latent HIV-1 reservoir is comprised of stably integrated and intact proviruses with limited to no viral transcription. It has been proposed that latent infection may be maintained by methylation of pro-viral DNA. Here, for the first time, we investigate the cytosine methylation of a replication competent provirus (AMBI-1) found in a T cell clone in a donor on antiretroviral therapy (ART). Methylation profiles of the AMBI-1 provirus were compared to other proviruses in the same donor and in samples from three other individuals on ART, including proviruses isolated from lymph node mononuclear cells (LNMCs) and peripheral blood mononuclear cells (PBMCs). We also evaluated the apparent methylation of cytosines outside of CpG (i.e., CpH) motifs. We found no evidence for methylation in AMBI-1 or any other provirus tested within the 5' LTR promoter. In contrast, CpG methylation was observed in the env-tat-rev overlapping reading frame. In addition, we found evidence for differential provirus methylation in cells isolated from LNMCs vs. PBMCs in some individuals, possibly from the expansion of infected cell clones. Finally, we determined that apparent low-level methylation of CpH cytosines is consistent with occasional bisulfite reaction failures. In conclusion, our data do not support the proposition that latent HIV infection is associated with methylation of the HIV 5' LTR promoter.

Short Communication: HIV-DRLink: A Tool for Reporting Linked HIV-1 Drug Resistance Mutations in Large Single-Genome Data Sets Using the Stanford HIV Database.

The prevalence of HIV-1 drug resistance is increasing worldwide and monitoring its emergence is important for the successful management of populations receiving combination antiretroviral therapy. It is likely that pre-existing drug resistance mutations linked on the same viral genomes are predictive of treatment failure. Because of the large number of sequences generated by ultrasensitive single-genome sequencing (uSGS) and other similar next-generation sequencing methods, it is difficult to assess each sequence individually for linked drug resistance mutations. Several software/programs exist to report the frequencies of individual mutations in large data sets, but they provide no information on linkage of resistance mutations. In this study, we report the HIV-DRLink program, a research tool that provides resistance mutation frequencies as well as their genetic linkage by parsing and summarizing the Sierra output from the Stanford HIV Database. The HIV-DRLink program should only be used on data sets generated by methods that eliminate artifacts due to polymerase chain reaction recombination, for example, standard single-genome sequencing or uSGS. HIV-DRLink is exclusively a research tool and is not intended to inform clinical decisions.

Linked dual-class HIV resistance mutations are associated with treatment failure.

We hypothesized that HIV-1 with dual-class but not single-class drug resistance mutations linked on the same viral genome, present in the virus population before initiation of antiretroviral therapy (ART), would be associated with failure of ART to suppress viremia. To test this hypothesis, we utilized an ultrasensitive single-genome sequencing assay that detects rare HIV-1 variants with linked drug resistance mutations (DRMs). A case (ART failure) control (nonfailure) study was designed to assess whether linkage of DRMs in pre-ART plasma samples was associated with treatment outcome in the nevirapine/tenofovir/emtricitabine arm of the AIDS Clinical Trials Group A5208/Optimal Combined Therapy After Nevirapine Exposure (OCTANE) Trial 1 among women who had received prior single-dose nevirapine. Ultrasensitive single-genome sequencing revealed a significant association between pre-ART HIV variants with DRMs to 2 drug classes linked on the same genome (dual class) and failure of combination ART with 3 drugs to suppress viremia. In contrast, linked, single-class DRMs were not associated with ART failure. We conclude that linked dual-class DRMs present before the initiation of ART are associated with ART failure, whereas linked single-class DRMs are not.

No evidence of HIV replication in children on antiretroviral therapy.

It remains controversial whether current antiretroviral therapy (ART) fully suppresses the cycles of HIV replication and viral evolution in vivo. If replication persists in sanctuary sites such as the lymph nodes, a high priority should be placed on improving ART regimes to target these sites. To investigate the question of ongoing viral replication on current ART regimens, we analyzed HIV populations in longitudinal samples from 10 HIV-1-infected children who initiated ART when viral diversity was low. Eight children started ART at less than ten months of age and showed suppression of plasma viremia for seven to nine years. Two children had uncontrolled viremia for fifteen and thirty months, respectively, before viremia suppression, and served as positive controls for HIV replication and evolution. These latter 2 children showed clear evidence of virus evolution, whereas multiple methods of analysis bore no evidence of virus evolution in any of the 8 children with viremia suppression on ART. Phylogenetic trees simulated with the recently reported evolutionary rate of HIV-1 on ART of 6 × 10-4 substitutions/site/month bore no resemblance to the observed data. Taken together, these data refute the concept that ongoing HIV replication is common with ART and is the major barrier to curing HIV-1 infection.

Ultrasensitive single-genome sequencing: accurate, targeted, next generation sequencing of HIV-1 RNA.

BACKGROUND: Although next generation sequencing (NGS) offers the potential for studying virus populations in unprecedented depth, PCR error, amplification bias and recombination during library construction have limited its use to population sequencing and measurements of unlinked allele frequencies. Here we report a method, termed ultrasensitive Single-Genome Sequencing (uSGS), for NGS library construction and analysis that eliminates PCR errors and recombinants, and generates single-genome sequences of the same quality as the "gold-standard" of HIV-1 single-genome sequencing assay but with more than 100-fold greater depth. RESULTS: Primer ID tagged cDNA was synthesized from mixtures of cloned BH10 wild-type and mutant HIV-1 transcripts containing ten drug resistance mutations. First, the resultant cDNA was divided and NGS libraries were generated in parallel using two methods: uSGS and a method applying long PCR primers to attach the NGS adaptors (LP-PCR-1). Second, cDNA was divided and NGS libraries were generated in parallel comparing 3 methods: uSGS and 2 methods adapted from more recent reports using variations of the long PCR primers to attach the adaptors (LP-PCR-2 and LP-PCR-3). Consistently, the uSGS method amplified a greater proportion of cDNAs, averaging 30% compared to 13% for LP-PCR-1, 21% for LP-PCR-2 and 14% for LP-PCR-3. Most importantly, when the uSGS sequences were binned according to their primer IDs, 94% of the bins did not contain PCR recombinant sequences versus only 55, 75 and 65% for LP-PCR-1, 2 and 3, respectively. Finally, when uSGS was applied to plasma samples from HIV-1 infected donors, both frequent and rare variants were detected in each sample and neighbor-joining trees revealed clusters of genomes driven by the linkage of these mutations, showing the lack of PCR recombinants in the datasets. CONCLUSIONS: The uSGS assay can be used for accurate detection of rare variants and for identifying linkage of rare alleles associated with HIV-1 drug resistance. In addition, the method allows accurate in-depth analyses of the complex genetic relationships of viral populations in vivo.

PAPNC, a novel method to calculate nucleotide diversity from large scale next generation sequencing data.

Estimating viral diversity in infected patients can provide insight into pathogen evolution and emergence of drug resistance. With the widespread adoption of deep sequencing, it is important to develop tools to accurately calculate population diversity from very large datasets. Current methods for estimating diversity that are based on multiple alignments are not practical to apply to such data. In this study, the authors report a novel method (Pairwise Alignment Positional Nucleotide Counting, PAPNC) for estimating population diversity from 454 sequence data. The diversity measurements determined using this method were comparable to those calculated by average pairwise difference (APD) of multiply aligned sequences using MEGA5. Diversities were estimated for 9 patient plasma HIV samples sequenced with Titanium 454 technology and by single-genome sequencing (SGS). Diversities calculated from deep sequencing using PAPNC ranged from 0.002 to 0.021 while APD measurements calculated from SGS data ranged proximately from 0.001 to 0.018, with the difference being attributable to PCR error (contributing background diversity of 0.0016 in a control sample). Comparison of APDs estimated from 100 sets of sequences drawn at random from 454 generated data and from corresponding SGS data showed very close correlation between the two methods with R(2) of 0.96, and differing on average by about 1% (after correction for PCR error). The authors have developed a novel method that is good for calculating genetic diversities for large scale datasets from next generation sequencing. It can be implemented easily as a function in available variation calling programs like SAMtools or haplotype reconstruction software for nucleotide genetic diversity calculation. A Perl script implementing this method is available upon request.

Low-frequency nevirapine (NVP)-resistant HIV-1 variants are not associated with failure of antiretroviral therapy in women without prior exposure to single-dose NVP.

BACKGROUND: Low-frequency nevirapine (NVP)-resistant variants have been associated with virologic failure (VF) of initial NVP-based combination antiretroviral therapy (cART) in women with prior exposure to single-dose NVP (sdNVP). We investigated whether a similar association exists in women without prior sdNVP exposure. METHODS: Pre-cART plasma was analyzed by allele-specific polymerase chain reaction to quantify NVP-resistant mutants in human immunodeficiency virus-infected African women without prior sdNVP who were starting first-line NVP-based cART in the OCTANE/A5208 trial 2. Associations between NVP-resistant mutants and VF or death were determined and compared with published results from women participating in the OCTANE/A5208 trial 1 who had taken sdNVP and initiated NVP-based cART. RESULTS: Pre-cART NVP-resistant variants were detected in 18% (39/219) of women without prior sdNVP exposure, compared to 45% (51/114) with prior sdNVP exposure (P < .001). Among women without prior sdNVP exposure, 8 of 39 (21%) with NVP-resistant variants experienced VF or death vs 31 of 180 (17%) without such variants (P = .65); this compares with 21 of 51 (41%) vs 9 of 63 (14%) among women with prior exposure (P = .001). CONCLUSIONS: The risk of VF on NVP-based cART from NVP-resistant variants differs between sdNVP-exposed and -unexposed women. This difference may be driven by drug-resistance mutations emerging after sdNVP exposure that are linked on the same viral genome. CLINICAL TRIALS REGISTRATION: NCT00089505.

Analysis of 454 sequencing error rate, error sources, and artifact recombination for detection of Low-frequency drug resistance mutations in HIV-1 DNA.

BACKGROUND: 454 sequencing technology is a promising approach for characterizing HIV-1 populations and for identifying low frequency mutations. The utility of 454 technology for determining allele frequencies and linkage associations in HIV infected individuals has not been extensively investigated. We evaluated the performance of 454 sequencing for characterizing HIV populations with defined allele frequencies. RESULTS: We constructed two HIV-1 RT clones. Clone A was a wild type sequence. Clone B was identical to clone A except it contained 13 introduced drug resistant mutations. The clones were mixed at ratios ranging from 1% to 50% and were amplified by standard PCR conditions and by PCR conditions aimed at reducing PCR-based recombination. The products were sequenced using 454 pyrosequencing. Sequence analysis from standard PCR amplification revealed that 14% of all sequencing reads from a sample with a 50:50 mixture of wild type and mutant DNA were recombinants. The majority of the recombinants were the result of a single crossover event which can happen during PCR when the DNA polymerase terminates synthesis prematurely. The incompletely extended template then competes for primer sites in subsequent rounds of PCR. Although less often, a spectrum of other distinct crossover patterns was also detected. In addition, we observed point mutation errors ranging from 0.01% to 1.0% per base as well as indel (insertion and deletion) errors ranging from 0.02% to nearly 50%. The point errors (single nucleotide substitution errors) were mainly introduced during PCR while indels were the result of pyrosequencing. We then used new PCR conditions designed to reduce PCR-based recombination. Using these new conditions, the frequency of recombination was reduced 27-fold. The new conditions had no effect on point mutation errors. We found that 454 pyrosequencing was capable of identifying minority HIV-1 mutations at frequencies down to 0.1% at some nucleotide positions. CONCLUSION: Standard PCR amplification results in a high frequency of PCR-introduced recombination precluding its use for linkage analysis of HIV populations using 454 pyrosequencing. We designed a new PCR protocol that resulted in a much lower recombination frequency and provided a powerful technique for linkage analysis and haplotype determination in HIV-1 populations. Our analyses of 454 sequencing results also demonstrated that at some specific HIV-1 drug resistant sites, mutations can reliably be detected at frequencies down to 0.1%.

Ultrasensitive allele-specific PCR reveals rare preexisting drug-resistant variants and a large replicating virus population in macaques infected with a simian immunodeficiency virus containing human immunodeficiency virus reverse transcriptase.

It has been proposed that most drug-resistant mutants, resulting from a single-nucleotide change, exist at low frequency in human immunodeficiency virus type 1 (HIV-1) and simian immunodeficiency virus (SIV) populations in vivo prior to the initiation of antiretroviral therapy (ART). To test this hypothesis and to investigate the emergence of resistant mutants with drug selection, we developed a new ultrasensitive allele-specific PCR (UsASP) assay, which can detect drug resistance mutations at a frequency of ≥0.001% of the virus population. We applied this assay to plasma samples obtained from macaques infected with an SIV variant containing HIV-1 reverse transcriptase (RT) (RT-simian-human immunodeficiency [SHIV](mne)), before and after they were exposed to a short course of efavirenz (EFV) monotherapy. We detected RT inhibitor (RTI) resistance mutations K65R and M184I but not K103N in 2 of 2 RT-SHIV-infected macaques prior to EFV exposure. After three doses over 4 days of EFV monotherapy, 103N mutations (AAC and AAT) rapidly emerged and increased in the population to levels of ∼20%, indicating that they were present prior to EFV exposure. The rapid increase of 103N mutations from <0.001% to 20% of the viral population indicates that the replicating virus population size in RT-SHIV-infected macaques must be 10(6) or more infected cells per replication cycle.

Short-course Combivir after single-dose nevirapine reduces but does not eliminate the emergence of nevirapine resistance in women.

BACKGROUND: In the Treatment Options Preservation Study (TOPS) trial, 4 or 7 days of Combivir (CBV; zidovudine/lamivudine) with maternal single-dose nevirapine (sdNVP) significantly reduced the emergence of NVP resistance as determined by virus population genotyping. To detect NVP resistance with greater sensitivity, we analysed TOPS samples by allele-specific real-time PCR (ASP). METHODS: In a random subset of women from each arm of the trial, plasma samples from before and 6 weeks after sdNVP were analysed using ASP at codons 103, 181, 184 and 190. RESULTS: Samples were analysed from 27 women in the sdNVP arm and 24 each in the CBV 4-day (sdNVP/CBV4) and 7-day (sdNVP/CBV7) arms. ASP detected NVP-resistant variants in week 6 samples from 70% of women in the sdNVP arm, 29% in the sdNVP/CBV4 arm and 33% in sdNVP/CBV7 arm (P<0.01 for sdNVP/CBV4 or sdNVP/CBV7 versus sdNVP; P=1.0 for sdNVP/CBV4 versus sdNVP/CBV7). Lamivudine resistance was detected by ASP in only 1 of 51 women who received CBV. CONCLUSIONS: Short-course CBV significantly reduced but did not eliminate the emergence of NVP resistance after sdNVP. NVP-resistant variants were detected in about one-third of women despite CBV treatment, but the duration of persistence and clinical impact of these variants in response to antiretroviral therapy is uncertain.

Role of low-frequency HIV-1 variants in failure of nevirapine-containing antiviral therapy in women previously exposed to single-dose nevirapine.

In the OCTANE/A5208 study of initial antiretroviral therapy (ART) in women exposed to single-dose nevirapine (sdNVP) ≥ 6 mo earlier, the primary endpoint (virological failure or death) was significantly more frequent in the NVP-containing treatment arm than in the lopinavir/ritonavir-containing treatment arm. Detection of NVP resistance in plasma virus at study entry by standard population genotype was strongly associated with the primary endpoint in the NVP arm, but two-thirds of endpoints occurred in women without NVP resistance. We hypothesized that low-frequency NVP-resistant mutants, missed by population genotype, explained excess failure in the NVP treatment arm. Plasma samples from 232 participants were analyzed by allele-specific PCR at study entry to quantify NVP-resistant mutants down to 0.1% for 103N and 190A and to 0.3% for 181C. Of 201 women without NVP resistance by population genotype, 70 (35%) had NVP-resistant mutants detected by allele-specific PCR. Among these 70 women, primary endpoints occurred in 12 (32%) of 38 women in the NVP arm vs. 3 (9%) of 32 in the lopinavir/ritonavir-containing arm (hazard ratio = 3.84). The occurrence of a primary endpoint in the NVP arm was significantly associated with the presence of K103N or Y181C NVP-resistant mutations at frequencies >1%. The risk for a study endpoint associated with NVP-resistant mutant levels did not decrease with time. Therefore, among women with prior exposure to sdNVP, low-frequency NVP-resistant mutants were associated with increased risk for failure of NVP-containing ART. The implications for choosing initial ART for sdNVP-exposed women are discussed.

Low frequency nonnucleoside reverse-transcriptase inhibitor-resistant variants contribute to failure of efavirenz-containing regimens in treatment- experienced patients.

BACKGROUND: The contribution of low frequency drug-resistant human immunodeficiency virus type 1 (HIV-1) variants to failure of antiretroviral therapy is not well defined in treatment-experienced patients. We sought to detect minor nonnucleoside reverse-transcriptase inhibitor (NNRTI)-resistant variants at the initiation of multidrug efavirenz-containing therapy in both NNRTI-naive and NNRTI-experienced patients and to determine their association with virologic response. METHODS: Plasma samples at entry and at time of virologic failure from patients enrolled in the AIDS Clinical Trials Group study 398 were analyzed by standard genotype, single-genome sequencing and allele-specific polymerase chain reaction (K103N and Y181C) to detect and quantify minor NNRTI-resistant variants. RESULTS: Minor populations of NNRTI-resistant variants that were missed by standard genotype were detected more often at study entry in NNRTI-experienced patients than NNRTI-naive patients by both single-genome sequencing (8 of 12 vs 3 of 15; P = .022) and allele-specific polymerase chain reaction (11% Y181C, 5 of 22 vs 3 of 72, respectively; P = .016). K103N variants at frequencies 11% were associated with inferior HIV-1 RNA response to efavirenz-containing therapy between entry and week 24 (change in HIV-1 RNA level, +0.5 vs -1.1 log(10) copies/mL; P < .001). CONCLUSIONS: Minor NNRTI-resistant variants were more prevalent in NNRTI-experienced patients and were associated with reduced virologic response to efavirenz-containing multidrug regimens.

Optimization of allele-specific PCR using patient-specific HIV consensus sequences for primer design.

Allele-specific PCR based on subtype consensus sequences is a powerful technique for detecting low frequency drug resistant mutants in HIV-1 infected patients. However, this approach can be limited by genetic variation in the region complementary to the primers, leading to variability in allele detection. The goals of this study were to quantify this effect and then to improve assay performance.

Suppression of viremia and evolution of human immunodeficiency virus type 1 drug resistance in a macaque model for antiretroviral therapy.

Antiretroviral therapy (ART) in human immunodeficiency virus type 1 (HIV-1)-infected patients does not clear the infection and can select for drug resistance over time. Not only is drug-resistant HIV-1 a concern for infected individuals on continual therapy, but it is an emerging problem in resource-limited settings where, in efforts to stem mother-to-child-transmission of HIV-1, transient nonnucleoside reverse transcriptase inhibitor (NNRTI) therapy given during labor can select for NNRTI resistance in both mother and child. Questions of HIV-1 persistence and drug resistance are highly amenable to exploration within animals models, where therapy manipulation is less constrained. We examined a pigtail macaque infection model responsive to anti-HIV-1 therapy to study the development of resistance. Pigtail macaques were infected with a pathogenic simian immunodeficiency virus encoding HIV-1 reverse transcriptase (RT-SHIV) to examine the impact of prior exposure to a NNRTI on subsequent ART comprised of a NNRTI and two nucleoside RT inhibitors. K103N resistance-conferring mutations in RT rapidly accumulated in 2/3 infected animals after NNRTI monotherapy and contributed to virologic failure during ART in 1/3 animals. By contrast, ART effectively suppressed RT-SHIV in 5/6 animals. These data indicate that suboptimal therapy facilitates HIV-1 drug resistance and suggest that this model can be used to investigate persisting viral reservoirs.

Blinded, multicenter comparison of methods to detect a drug-resistant mutant of human immunodeficiency virus type 1 at low frequency.

We determined the abilities of 10 technologies to detect and quantify a common drug-resistant mutant of human immunodeficiency virus type 1 (lysine to asparagine at codon 103 of the reverse transcriptase) using a blinded test panel containing mutant-wild-type mixtures ranging from 0.01% to 100% mutant. Two technologies, allele-specific reverse transcriptase PCR and a Ty1HRT yeast system, could quantify the mutant down to 0.1 to 0.4%. These technologies should help define the impact of low-frequency drug-resistant mutants on response to antiretroviral therapy.